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ABSTRACT 

We present the analysis of the radio jet evolution of the radio galaxy 3C 120 during a period of 
prolonged 7-ray activity detected by the Fermi satellite between December 2012 and October 2014. We 
find a clear connection between the 7-ray and radio emission, such that every period of 7-ray activity 
is accompanied by the flaring of the mm-VLBI core and subsequent ejection of a new superluminal 
component. However, not all ejections of components are associated with 7-ray events detectable 
by Fermi. Clear 7 — ray detections are obtained only when components are moving in a direction 
closer to our line of sight. This suggests that the observed 7-ray emission depends not only on the 
interaction of moving components with the mm-VLBI core, but also on their orientation with respect 
to the observer. Timing of the 7-ray detections and ejection of superluminal components locate the 
7-ray production to within ^0.13 pc from the mm-VLBI core, which was previously estimated to lie 
about 0.24 pc from the central black hole. This corresponds to about twice the estimated extension 
of the broad line region, limiting the external photon field and therefore suggesting synchrotron self 
Compton as the most probable mechanism for the production of the 7-ray emission. Alternatively, 
the interaction of components with the jet sheath can provide the necessary photon field to produced 
the observed 7-rays by Compton scattering. 

Subject headings: galaxies: active — galaxies: jet — galaxies: individual (3C120) — radio continuum: 
galaxies 


1. INTRODUCTION 

In the unified scheme of active galactic nuclei (AGN), 
Fanaroff-Riley radio galaxies of type I (FRI) and 11 
(FRH) are considered the parent population of BL Lacs 
and flat spectrum radio quasars (FSRQ), respectively. 
FRI and FRH radio galaxies belong to the misaligned 
AGN class, as they are oriented at larger viewing angles 
than blazars (the most luminous and variable BL Lac 
objects and flat-spectrum radio quasars). Relativistic 
beaming amplifies the emission of jets pointing toward 
the observer, making blazars the brightest objects in the 
extragalactic sky also in the 7-ray band. In fact, among 
the more than one thousand extragalactic sources de¬ 
tected by the Large Area Telescope (LAT) onboard the 
Fermi Gamma-Ray Space Telescope in two years, only 
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3% of them are not associated to blazar objects. FRI and 
FRH radio galaxies fall inside this 3%, wit h the predom¬ 


inance of nearby radio galaxies of type I (Grandi et al. 


2012a). FRH radio galaxies with GeV emission are rare 
being 3G 111 the only FRH so far with a confirmed Fermi 
counterpart. 3G 111 is also the first FRH radio galaxy 
where a 7-ray flare has been associated with t he ejec- 
tion of a new bright knot from the radio core ( Grandi] 


et al. 2012bI. The simultaneity of the observed hare in 
millirneter, op tical, X-ray and 7-ray bands lead |Grandi| 


et al. (2012b) to claim that the GeV dissipation region 
is located at a distance of about 0.3 pc from the central 
black hole. 

The rapid 7-ray variability observed in FRI ra¬ 
dio galaxies , with time scales of months in the cas e of 


NGG 1275 (Abdo et al. 2010 Kataoka et al. 2010), or 


even days for the case of the I'eV vari abilit-y in M87 
( Aharonian et al.|[2006 Harris et al.|2011 ), suggests that 
the 7-ray emission m these sources originates also in a 
very compact jet region. Results fro m over ten years of 


multi wavelength observations of M87 (Abramowski et al. 
2012[) suggest that the very high energy hares may take 


place in the core (for the case of the 2008 and 2011 TeV 
flares), or in the HST-1 complex, about 0.8" downstream 
the jet (in the cas e of the 2005 TeV flare). More recently, 


Hada et al. (|2014) found that the last TeV flare occurred 
in iVI87 in 2012 originates in the jet base, within 0.03 pc 
from the black hole, but no correlation with the MeV- 
GeV light curve obtained by Fermi has been found. 

Recently, the Large Area Telescope detected an un¬ 
precedented 7-ray flare from the FRI radio galaxy 3G 120 
on September 24th, 2014 when the source reached a daily 
flux (E>100 MeV) of 1.0±0.3xl0~^ photon cm“^ s“^ 


as reported in Tanaka et al. (2014). This recent flare 
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appears to be associated with a higher state of 7 -ray ac¬ 
tivity in the source. In fact, since December 2012 the 
radio galaxy 3C 120 registered a series of 7 -ray events, 
indicating a flaring activity that lasted until at least our 
last data analyzed, in 2014 October 4. 

The radio galaxy 3C 120 (z=0.033) presents an FRI 
morphology, but it also has a blazar-like radio jet, show- 
ing multiple superlumina l comp onents at parsec scales 


(Gomez et al. 2000 2001 20111, as well as at distances 


as large as 150 pc in projection from the core (Walker 


et al. 2001). This radio galaxy reveals also X-ray prop- 
erties similar to Seyfert galaxies, with the X-ray spectra l 
slope increasing with intensity (Maraschi et al. 1991), 


and a prominent iron emission line at a photon energy 
of 6.4 KeV. This implies that most of the X-ray emission 
comes from, or near the accretion disk, rather than in the 
jet. In addition, the observed strong correlation between 
dips in the X-ray emission and the ejec tion of new su¬ 
perluminal components in th e radio jet (Marscher et al. 


2002 Chatterjee et al. 


2009) reveals a clear connection 


between the accretion disk and the radio jet. 

In this paper we present the first association of 7 -ray 
emission and the ejection of new superluminal compo¬ 
nents in a FRI radio galaxy, 3C 120. This resem bles the 
recent findin gs on the FRII radio galaxy 3C 111 (Grand! 
et al.||2012b). We present the radio data set analyzed in 
this study, as well as methods used to reduce radio data 
in Section 2; in Section 3 we present the analysis and 
results of the Fermi-LKI data; in Section 4 we study the 
radio emission at 15 and 43 GHz in the parsec scale jet; 
Section 5 presents the connection between the 7 -ray and 
radio emission, and in Section 6 we discuss our findings. 

The cosmological values adopted from the Planck’s re¬ 
sults ( Planck Collaboration et al.| 2014) are 11^= 0.3, 
r2A= 0.7, and Hq = 68 km s“^Mpc"^. With these values, 
at the redshift of 3C 120 (z=0.033) 1 mas corresponds to 
a linear distance 0.67 pc, and a proper motion of 1 mas 
yr“^ corresponds to an apparent speed of 2 . 21 c. 

2. RADIO DATA ANALYSIS 

To study the structure of the radio jet in 3C 120 we 
have collected data from two of the most extended VLBA 
monitoring programs: the MOJAVEp^ and the VLBA- 
BU-BLAZAR programj^ This radiodataset consists of 
46 epochs of VLBA data at 15 GHz taken from the MO¬ 
JAVE survey, covering the observing period from June 
2008 to August 2013, and 21 epochs of VLBA data at 43 
GHz from the VLBA-BU-BLAZAR program, covering 
the period from January 2012 to May 2014. 

The reduction of the VLBA 43 GHz data has been 
performed using a comb ination of AIPS and Difmap 
packages, as described in Jorstad et al. (2005). VLBA 


data at 15 GHz have been calibrated by t he MOJAVE 
team, f ollowing the procedure described in [Lister et al.| 
(2009a). Eor comparison across epochs all the images 
have been convolved with a mean beam of 0.3x0.15 mas 
and 1.2X0.5 mas for the VLBA-BU-BLAZAR and MO¬ 
JAVE programs, respectively. 

To determine the structural changes in the radio jet we 
have modeled the radio emission through fitting of the 
visibilities to circular Gaussian components using Difmap 


(Shepherd 1997). Fitted values for each component are 
the hux density, separation and position angle from the 
core, and size. These are tabulated in two tables at the 
end of the paper - Tables and [i] for the 15 and 43 GHz 
data, respectively. 

3. FERMI-LAT DATA ANALYSIS 

The LAT data collected during 72 months of opera¬ 
tion (from 2008 August 4 to 2014 August 4were an¬ 
alyzed using the Fermi-LKT ScienceTools software (ver¬ 
sion v9r32p5) and the P7REP_SOURCE_V15 set of in¬ 
strument response functions (Ackermann et al.||2012||^^ 
The time intervals when the rocking angle of the LAL 
was greater than 52° were rejected and a cut to select a 
maximum zenith angle of 100 ° of the events was applied 
to exclude 7 -rays originating from cosmic ray interactions 
with the Earth’s atmosphere. 

The detection significance of a source is provided by the 
TS = 2[logL(source)-logL(no source)], where L(source) 
is the maximum likelihood value for a model with an 
additional source at a specified location and L(no source) 
is the maximum l ikelihood value for a model without the 
additional source(Mattox et al. 1996). When the TS is 


http://www.physics.purdue.edu/MOJAVE/ 
http://www.bu.edu/blazars/research.html 


less than 10 or the ratio of the hux uncertainty to the 
flux is more than 0.5, a 2cr upper limit of the flux is 
provided. Depending on the TS value, the upper limits 
are calculated using the profile (TS> 1) or the Bayesian 
(TS< 1) metho d as described in t he second LAT catalog 
(2FGL catalog; Nolan et al.|2012 ). All errors reported in 
the figures or quoted in the text are la statistical errors. 

The estimated systematic errors on the flux, 10% at 100 
MeV, decreasing to 5% at 560 MeV, and increasing to 
10% at 10 GeV, refer to uncertainties on the effective 
area of the instrumenlF^ 

3.1. The 72-month Average Spectrum 

We performed both binned and unbinned likeli¬ 
hood analysis following the standard LAT data anal¬ 
ysis procedures, obtaining consistent results. Here 
we present results from binne d analysis to be 
consi stent with published data (Ackermann et al. 

2011). We accumulated 72 months of data to ob- 
tam the average 7 -ray spectral properties of the 
source. The adopted model included all 2EGL 

sources within 15° of 3C 120 (R.A.(J2000)=68.2962313°, 
Dec(J2000)=5.3543389°). The studied source was mod¬ 
eled with a power law. All spectral parameters of the 
sources more than 10° from the center of the Region of 
Interest (Rol) were fixed to the 2EGL values. The Galac¬ 
tic diffuse emission was modeled using the standard dif¬ 
fuse emission model glLiem_v05-revl.fit while isotropic 
7 -ray emission and the residual cosmic ray contamina¬ 
tion in the instrument were modeled using the template 
iso_source-v05. txt. 

3G 120 was detected in the 100 MeV-100 GeV band 
with a TS value of 107 (^ IOct). The source is steep 
( photon index r=2.7±0.1) and weak with a flux of 

Mission Elapsed Time (MET) Start Time=239557417:MET 
End Time = 428803203 

Science Tools and instrument response functions 
are available from the Fermi Science Support Center: 
http://fermi.gsfc.nasa.gov/ssc/data/analysis. 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats_pass7.html. 
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Fig. 1.— Fermi-LAT light curves of 3C 120 covering 72-month of survey (from 2008 August 4 to 2014 August 6 ) obtained considering 
two different energy bands: 100 MeV-100 GeV {left panel) and 500 MeV-100 GeV {right panel). A bin width of 3 months is adopted. 
A 2 cr upper limit flux is shown (arrow) when the source is not detected (TS<10). For comparison the light curves of the nearby FSRQ 
(PKS 0423-1-051) at 1.6 degrees from 3C 120 are also shown in the same energy ranges. There is not significant overlapping, as 3C 120 and 
PKS 0423-h051 appear to be active in different time intervals. 


F>ioo MeV = 2.5 ± 0.4 X 10 ^ photon cm 




agreement with the result of Abdo et al. (2010). As 
noted by the same authors, a h'iat Spectrum tjuasar, 
PKS 0423+05, is located at only ~ 1.6 degrees from 
the radio galaxy. In our analysis, the blazar is soft 
(r=2.6+0.1) and slightly fainter (F>ioo MeV = 1-77 ± 
0.36) X 10“® photon cm“^ s“^) than 3C 120. 

3.2. Light Curves 

Gamma-ray light curves were produced by dividing the 
analyzed time interval in temporal segments and repeat¬ 
ing the likelihood analysis with only the normalizations 
of the sources within 10° free to vary. The spectral slopes 
of all the sources in the Rol were kept fixed to the best 
fit values of the 72-month likelihood analysis. 

At first we produced a light curve from 2008 August 4 
to 2014 August 4 with a bin size of 3 months in the 100 
MeV-100 GeV energy band. The light curve shown in 
Figured (left upper pan el) indicates that 3 C 120 was ac¬ 
tive starting from 2012 . Abdo et al. (2010) reported also 
a detection between November 2008 and h'ebruary 2009. 
Our new analysis, which was performed with up-dated 
background files and new instrument response functions, 
provides only a flux upper limit for the same time in¬ 
terval (second bin of the light curve). Although we can 
not claim a detection, the calculated TS value is however 
high (TS=9.2, corresponding to ^ 3cr) and only slightly 
below the usually adopted detection threshold (TS=10). 

Gonsidering that PKS 0423+051 is at only 1.6° from 
3G120, that both sources have steep spectra and, that 
the LAT PSF is very broad at MeV energies, we de¬ 
cided to explore the 3-month light curve of the quasar to 
check for possible confusion effects. As shown in Figured 
PKS 0423+051 and 3C 120 flare in different periods and 
exhibit light curves with different pattern, suggesting a 


negligible (mutual) contamination. 

As further test, we also performed a variability study 
considering only photons with energy > 500 MeV. Above 
this energy, the 68 % containment angle (i.e., the radius 
of the circle containing 68 % of the PSF) is indeed compa¬ 
rable or smaller than the separation of the two sources. 
The 500 MeV—100 GeV light curves (Figure right 
panel) are similar to those obtained taking also into ac¬ 
count softer photons, supporting the previous conclusion 
on negligible confusion effects. 

Finally, we reduced the integration time interval of 
each bin to 15 days to better constrain the time of the 
7 -ray detections, shown in Table In this case the 
analysis has been extended until October 2014 to con- 
firm the flare at the end of September 2014 reported in 


Tanaka et al. (2014). Our analysis yields a 7 —ray flux of 
2.52 +0.86 X 10 ® photon cm 


-2 s-i in MJD 56924-56939 


for the September 2014 event (see Table [^. 


TABLE 1 

7-RAY DETECTIONS OF 3C 120 IN THE ENERGY RANGE 0.5-100 
GeV with a bin size of 15 daysH 

^Notes.We report the integration time in MJD and date, the flux 
with the corresponding error in 10“® photon cm~^ and the TS 
value associated to each 7 -ray detection. 


MJD 

Date 

Flux 

Err 

TS 

56264-56279 

03/12/2012-18/12/2012 

1.41 

0.60 

10.6 

56549-56564 

14/09/2013-29/09/2013 

1.65 

0.54 

21.6 

56564-56579 

29/09/2013-14/10/2013 

1.31 

0.53 

13.3 

56774-56789 

27/04/2014-12/05/2014 

1.15 

0.52 

10.4 

56819-56834 

11/06/2014-26/06/2014 

1.47 

0.59 

13.8 

56924-56939 

24/09/2014-09/10/2014 

2.52 

0.86 

18.4 
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4. THE PARSEC SCALE JET AT 15 AND 43 GHZ 
4.1. VLB A data at 15 GHz 

The study of the jet evolution at 15 GHz has been 
performed on the series of 46 VLBA images obtained by 
the MOJAVE monitoring progr^am - a subset of these 
images is displayed in Figure where contours repre¬ 
sent the total intensity with model-fit components (red 
circles) overlaid. 

We detected in total 25 components, apart from the 
core, although some of them characterize more probably 
the underlying flux density than knots that have been 
ejected from the core and move along the jet. The radio 
core is usually defined as the bright, compact feature 
at the upstream end of the jet, which may correspond 
to a recollimation shock (i.e., Gomez et al. 1997) at 
millimeter wavelengths and to the optically thin-thick 
transition at centimeter wavelengths. 

Components E-F and C are robustly identified mov¬ 
ing and standing features, respectively, and the rest of 
model-fit components are required by the data, but can¬ 
not be cross-identified across the epochs and may rep¬ 
resents e.g., emission from the underlying jet flow. The 
core, identified with component CO, is considered sta¬ 
tionary across epochs. Plots of the separation from the 
core and flux density evolution of the fitted components 
are shown in Figs. and respectively. 

We fit the trajectories with respect to the core for all 
the superluminal knots that can be followed for a signif¬ 
icant number of epochs, as plotted in Fig. In order 
to have a better determination of the time of ejection 
of each component, namely the time when a new knot 
crosses the radio core, we use linear fits for the compo¬ 
nent separation versus time. Note that in some cases we 
have found evidence for components merg ing, splitting, 


or a clear a cceleration in their motion (e.g., Homan et al. 


2009 20151. In those cases we have considered only the 


initial epochs with a clear linear fit, as we are mainly 
interested in determining the time of ejection of each 
component. The time of ejection, angular and apparent 
velocities of the components are tabulated in Table 
We find that superluminal components move with ap- 


parent velocities between 5 and 6 c, in ag 

reement with 

previous findings ( 

Gomez et al. 12001 

2008 

2011 

Jorstad 

et al. 

2005). This 

agrees also with recent MUJ 

AVE re- 


sults (Lister et al. 2013) where the 3C 120 kinematics, 
together with the radio galaxy 3C 111, departs from the 
others in the sample. These two radio galaxies seem to 
be the only ones displaying clear superluminal motions 
and with apparent speeds that do not commonly change 
with distance from the core. 

From June 2008 to the beginning of 2012 a stationary 
feature. Cl, is found at a distance that shifts between 
^ 0.4 - 0.7 mas from the c ore. We identify this com po- 
nent with that reported by Leon-Tavares et al. (2010) at 
0.72±0.25 mas from the core, which is assumed by the 
authors to be related with optical flares when new com¬ 
ponents cross its position. Starting in 2012 the innermost 
1 mas region of the jet changes, so that two stationary 
components can be found, Cl and C2. An extra com¬ 
ponent, labeled S9, is also required to fit some extended 
and weak flux density at ^1.5 mas from the core. 

As can be seen in Figs.j^andffl the presence of two sta¬ 
tionary components (Cl and C2) within the innermost 


1 mas is associated with an increase in flux density and 
a more extended structure of the core. Component Cl 
is observed to progressively increase its flux density be¬ 
tween mid 2012 and the beginning of 2013, when both, 
the core and Cl are in a high flux density state. Three 
months later component C2 also shows an increase in flux 
density. These changes in the innermost structure of the 
jet are associated with the ejection of a new component, 
dll , revealed by the 43 GHz VLBA data (see section 


4.2), which extends our study of the jet until May 2014. 


wo other radio flares in the core have been also ob¬ 
served, one at the end of 2008 (2008.76±0.10 yr) and a 
second one at the beginning of 2011 (2011.15±0.05 yr). 
The flare in 2011 is associated with the ejection of a 
weak radio component from the core, E9, for which we 
estimate the time of ejection in 2011.23±0.04. 


TABLE 2 

Time of ejection, proper motions and apparent velocities 
OF MODELFIT COMPONENTS AT 15 GHz. 


Name 

T • 

(year) 

(mas/yr) 

^app 

EO 

2007.29 ± 0.06 

2.81 ± 0.05 

6.21 ± 0.11 

El 

2008.01 ± 0.02 

2.76 ± 0.05 

6.10 ± 0.11 

E4 

2008.82 ± 0.04 

2.35 ± 0.05 

5.19 ± 0.11 

E5 

2009.42 ± 0.02 

2.60 ± 0.06 

5.75 ± 0.13 

E6 

2009.85 ± 0.03 

2.56 ± 0.05 

5.66 ± 0.11 

E8 

2010.45 ± 0.02 

2.20 ± 0.03 

4.86 ± 0.07 

E9 

2011.23 ± 0.04 

2.32 ± 0.09 

5.12 ± 0.19 


The radio core flare at the end of 2008 is instead as¬ 
sociated with the ejection of a bright superluminal knot, 
E4, in 2008.82±0.04 yr. The light curves of the different 
components shown in Fig. reveal the unusually large 
flare experienced by E4. By mid 2009 it has doubled its 
flux density, reaching a peak of ~ 1.5 Jy and becoming 
brighter than the core itself. During the high state of flux 
density of E4 another component, named E4a, appears 
very close to it, suggesting an increase in extension of 
component E4 during the flare, as can be seen in Figs. 

and 1^ In this case we consider more appropriate to 
take into account a flux-density-weighted distance be¬ 
tween that of E4 and E4a to estimate the distance from 
the core of the new component E4 during its high state 
of flux density. The same method is used for compo¬ 
nents El and Ela, as we considered Ela an extension of 
El soon after this component is ejected from the core at 
the beginning of 2008. A similar splitting of components 
was also observed previousl y in 43 GHz VLBA images of 
3C 120 dComez et aL|[MI| ) 

The extended emission structure of components El and 
E4 (associ ated with compone nts Ela and E4a) is consis¬ 
tent with Aloy et al. (2003) relativistic hydrodynamic 
simulations, where the passage of a new perturbation 
from a series of recollimation shocks results in extended 
emitting regions due to light-travel time delays between 
the front and the back of the perturbation. 

Figure shows also the light-curve obtained from 
adding the flux densities of CO and Cl. The station¬ 
ary component Cl is usually located at a distance from 
the core of the order or smaller than the observing beam, 
therefore in many epochs it is difficult to disentangle its 
flux density from that of the core. For instance, the com- 
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Fig. 2. — Sequence of total intensity 15 GHz VLBA images from the MOJAVE monitoring program of 3C 120 with a common restoring 
beam of 1.2x0.5 mas at 0°. The separation among images is proportional to the time elapsed between observing epochs. Contours are 
traced at 0.0015, 0.004, 0.009, 0.02, 0.05, 0.1, 0.3, 0.6, 1.2 Jy. Red circles represent model-fit components. 


bination of the CO+Cl flux density reveals a high state 
leading to the ejection of component E 8 in mid 2010. It 
is also particularly remarkable the increase in flux den¬ 
sity of the CO-l-Cl complex in mid 2013 leading to the 
ejection of a new component, dll, seen at 43 GHz (see 
below). 

4.2. VLBA data at 43 GHz 

We modeled the 43 GHz jet in the same manner as 
for the 15 GHz data, finding a total of 16 components, 
although as in the case of the 15 GHz data, some of them 
are most probably related to the underlying continuum 


flow and it is difficult to follow them along the jet. In 
each epoch we identify the core with cO and we consider 
it a stationary feature. 

Figure shows a sequence of VLBA 43 GHz images 
in total intensity, where red circles represent model-fit 
components. The displayed images cover the observing 
period from October 2012 to May 2014, when 3C 120 
was detected at 7 -ray frequencies. 

As observed in the 15 GHz data, starting in mid 2012 
we find an increase in the flux density and extension of 
the core, requiring up to three stationary components. 
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2008 2009 2010 2011 2012 2013 



Fig. 3. — Distance from the core vs. time for the 15 GHz model-fit 
components with linear fits overlaid. Downward black arrows mark 
the time of ejection of each component. Grey vertical lines indicate 
7 -ray detections in the 15 days-bin (dark grey) and 3 months-bin 
(light grey) light curves in the energy band 500 MeV-1 00 GeV. The 
orange verti cal line marks the optical flare reported in |Kollatschny | 
|et al. 1 1 |2014[ | (See text). 


2008 2009 2010 2011 2012 2013 



Fig. 4.— Light curves of the 15 GHz model-fit components, in¬ 
cluding the added total flux density of the core and component G1 
(COH-Cl). Arrows and vertical lines are as indicated in Fig. [3| 



dlO, d9, and d^ to model the innermost 0.5 mas struc- 
ture (see Figs. and [^. Because of the different angular 
resolution and opacity with respect to that obtained at 
15 GHz, it is difficult to establish a one-to-one connection 
among the components at 43 and 15 GHz for the inner- 


Fig. 5.— Sequence of total intensity 43 GHz VLB A images from 
the VLBA-BU-BLAZAR monitoring program of 3C 120 with a 
common restoring beam of 0.3x0.15 mas at 0°. Red circles rep¬ 
resent modelfits components. Gontours are traced at 0.16, 0.35, 
0.77, 1 .7, 3.77, 8.33, 18.41, 40.71, 90% of the peak intensity, 1.82 
Jy/beam. 
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2012 2013 2014 



Fig. 6. — Distance from the core vs. time for the 43 GHz model-fit 
components with linear fits overlaid. Downward black arrows mark 
the time of ejection of each component. Grey vertical lines indicate 
7 -ray detections in the 15 days-bin (dark grey) and 3 months-bin 
(light grey) light curves in the energy band 500 MeV-100 GeV. 


2012 2013 2014 



the perturbation crosses the last stationary component, 
d 8 , we can clearly discern a new knot, dll, that emerges 
from the first 0.5 mas in 2013.4 yr (see Fig. [^. From the 
progressive flaring of stationary features close to the core 
we can infer the extension of the crossing emitting region, 
obtaining ~0.35 mas. This is significantly larger that the 
FWHM obtained from the model-fit of dll (see Table®, 
which suggests that dll is in fact part of a more extended 
region, resembling the results obtained from r elativistic 
hydrodynamical simulations (Aloy et al. 2003). Fitting 
of the separation versus distance for dll yields a proper 
motion of 1.91±0.09 mas/yr, which corresponds to an 
apparent velocity of 4.22±0.22c. The estimated time of 
ejection, that is, when component dll crossed the radio 
core at 43 GHz, is 2013.03±0.03. 

A similar situation takes place also in the second half 
of 2013, when the core, together with components dlO 
and d9, starts to bright and a new component, dl2, 
appears at ~ 0.45 mas. We note also that another 
component, dl 2 a, appears very close in time and po¬ 
sition to component dl 2 , although it displays a signif¬ 
icantly slower proper motion. Component dl2 moves 
at 2.1±0.2 mas/yr (4.7±0.3 c), while component dl2a 
moves at 1.2±0.2 mas/yr (2.6±0.5 c). We note, however, 
that the estimated time of ejection for dl 2 and dl 2 a, 
2013.67±0.02 yr and 2013.64±0.07 yr, respectively, is the 
same within the uncertainties. Both components there¬ 
fore originated simultaneously in the millimeter VLBI 
core, although they propagate at quite different veloci¬ 
ties afterwards. 

Time delays stretch the shocked emission in the ob¬ 
server’s frame, so that with the necessary angular res¬ 
olution multiple sub-components associated with a sin¬ 
gle shock could be distinguished in the jet, but in this 
case they would have similar apparent velocities, in con¬ 
trast with what it is observed for components dl 2 and 
dl2a. Trailing components have a smaller velocity than 
the leading perturbation, b ut they are release d on the 
wake of main perturbation (Agudo et al. [ 2001 ), instead 
of being ejected from the core, as occurs for component 
dl2a. On the other hand, relativistic hydrodynamic sim¬ 
ulations show that a single perturbation in the jet inlet 
leads to the formation of a forward and reverse shock 


(I Gomez et al. |1997t|Aloy et al.|2003||Mimica et al.|2009|). 

i'herefore, the fact that componeiits dl 2 and dl 2 a are 

ejected from the core at the same time but with differ¬ 
ent velocities suggests that they may correspond to the 
forward and reverse shock of a perturbation. 


Fig. 7.— Light curves of the 43 GHz model-fit components. Ar¬ 
rows and vertical lines are as indicated in Fig. 

most jet region. We tentatively identify components s9b 
and s9c with substructures of component S9 seen at 15 
GHz. Furthermore, during the enhanced activity of the 
source starting in mid 2012 it is not clear whether dlO, 
d9, and d 8 at 43 GHz - and Gl, G2 at 15 GHz - corre¬ 
spond to actual physical structures in the jet, like recol- 
limation shocks, or they trace the underlying emission of 
the jet. 

The radio core persists in a flaring state until January 
2013. After this, the peak in flux density moves from 
the core along the jet, crossing progressively the three 
stationary features close to it (dlO, d9 and d 8 ). When 


5. CONNEGTION BETWEEN 7 -RAY AND RADIO 
EMISSION 

Our an^sis of the Fermi-LAT data (see Figs. 
and Tab. [H shows a prolonged 7 -ray activity in the 
radio galaxy 3C 120 between December 2012 and 
September-October 2014, when the source reaches a flux 
ofF> 5 oo MeV = 2.5 ± 0.8 X 10“® photon cm“^ s“^, about 
an order of magnitude larger than in previous detections. 
Three clear periods of 7 -ray activity are found at the end 
of 2012 (56264-56279 MJD), September-October 2013 
(56549-56579 MJD), and May-October 2014 (56774- 
56939 MJD). 

The 7 -ray activity in December 2012 is accompanied 
by an increase in the radio core flux density, as well as 
in the innermost stationary components, at both 15 and 
































43 GHz (see Figs. Handl^, leading to the ejection of 
component dll in 2013.13^0.03 (see also Fig.I^. From 
the estimated proper motion of dll we infer that in the 
last epoch of the 15 GHz data this component is still 
in the innermost region of the jet as imaged at 15 GHz, 
crossing the stationary component C2. 

After this first 7 -ray detection the source remains in 
a quiescent state until August 2013, when another 7 -ray 
event covering the period from August 2013 to Octo¬ 
ber 2013 is detected. The 15-days binning analysis (see 
Tab.fH constrains the main flaring activity to September- 
OctoMr 2013 (56549-56579 MJD), coincident with the 
flaring of the radio core at 43 GHz and the ejection of 
the dl2-dl2a pair of components (see Figs. to [^. 

There is also indication of the beginning of a new flar¬ 
ing activity in the 43 GHz radio core associated with the 
last 7 -ray flaring activity starting in May 2014, but fur¬ 
ther VLBI images are required to confirm the ejection of 
a new component in this event. 

We therefore can conclude that there is a clear asso¬ 
ciation between the 7 -ray and radio emission in 3C 120, 
such that every period of 7 -ray activity is accompanied 
by flaring of the VLBI radio core and subsequent ejection 
of a new super luminal component in the jet. However, 
not all ejections of superluminal components are related 
to enhanced 7 -ray emission, detectable by Fermi-LAT, as 
occurred for components E4, E5, E 6 , E 8 , and E9 . The 
case of com ponent E4 is particularly interesting. |Abdo| 
et al. ( |2010 ) report a 7 -ray detection in December 2008, 
coincident with an increase of the 15 GHz radio core flux, 
leading to the ejection of E4 in 2008.82±0.04. This is also 
coincident with an optical fla r e in D ecember 2008, as re¬ 
ported by Kollatschny et al. (2014). We also note that 
E4 reached a peak hux of ~ 1.5 Jy, significantly larger 
than any other component seen in our analysis. Despite 
this intense activity in the optical and radio core, our 
analysis of the Fermi data during this period does not 
provide a clear detection in 7 -rays, although the calcu¬ 
lated TS (corresponding to ~ 3cr) is only slightly smaller 
than the usually adopted detection threshold. 

6. DISCUSSION 

6.1. Motion of components along the jet 

The debate on the origin and location of the 7 -ray 
emission in blazars has gained added interest since the 
launch of the Fermi satellite. Much of the current dis¬ 
cussion lies in whether 7 -rays are produced upstream of 
the mm-VLBI core, a s suggested by some 7 -ra y and ra¬ 
dio correlations (e.g. Rani et al. 2013a 20 14|) and the 
observed 7 -ray spectral break at few Ge v^(|Abdo et al. 
'20091 Finke fc Dermerl 2010[ Poutanen fc Stern '2010‘ 


Tanaka et al. 2011| Rani et al. 2013bp or downstream 

as suggested by coincidence of 7 -ray flares with eithei 


7 -ray 


either 


the appearance of new superluminal components (e.g. 
Jorstad et al.||2010 2013) or the passage of moving com¬ 


ponents through a st ationary jet feature (Schinzel et al. 


As observed in several blazars I 

Jorstad et al. 

2010 

2013 Ramakrishnan et al. 2014 

), the interaction be 

tween 


traveling features and the stationary radio core appears 
to be a necessary condition for the production of 7 -ray 
photons in 3G 120, but it is clearly not enough. There¬ 
fore, to understand the 7 -ray emission in 3G 120, and 


more generally in AGN, it is necessary to address the 
question of what physical changes in the jet can produce 
7 -ray emission. 

We note that the beginning of the 7 -ray activity in 
3G 120 occurs after a sustained period of low activity 
in the jet. As indicated in Table during the time 
period analyzed new components are seen in the jet of 
3G 120 roughly every 8 months; however, no new com¬ 
ponents are detected in the jet between the ejection of 
E9 in 2011.23±0.04 and component dll in 2013.03±0.03, 
implying a lack of activity in the jet for almost 2 years. 
Note that after the ejection of dll the source resumes 
its activity with the ejection of component dl 2 , again 
roughly 8 months later. 

Analysis of the components’ proper motions (see Ta¬ 
ble reveals a clear pattern of decreasing apparent ve¬ 
locities, from the 6.21±0.11c of EO to 5.12±0.19c of E9. 
Furthermore, when the core resumes the ejection of com¬ 
ponents after a 2 -year inactivity, the apparent veloci¬ 
ties measured for dll and dl 2 have further decreased 
to 4.22±0.22c and 4.70±0.31c, respectively. Gomponent 
dl2a shows an even smaller velocity of 2.75±0.45c, but 
we believe this component differs for the other compo¬ 
nents seen previously in 3C 120 in which it is probably 
associated with a reserve shock. 

Note that model-fits at 43 GHz and 15 GHz trace the 
motion of components continuously throughout the jet. 
We also find no evidence for acceleration (only compo¬ 
nent E 6 shows a clear acceleration at ~8 mas from the 
core), hence we consider that the proper motions mea¬ 
sured at 43 GHz provide a good estimation for the ex¬ 
pected values should these components be detected later 
on at 15 GHz, reassuring our finding for a progressive de¬ 
crease in the apparent velocity of the components from 
^ 6 .2c to 4.2c in a time span of approximately 6.4 years, 
from 2007.3 to 2013.7. This progressive change in the 
apparent velocity of components could be due to a change 
in the velocity and/or orientation of the components. 
Gonsidering that the 6.4 yr time span measured roughly 
agrees with the 12.3±0.3 yr full period determined by 
the precessing jet model of Caproni & Abraham (2004), 
we favor the case in which it is produced by a change 
in the orientation of the components with respect to the 
observer. 

Figure |8 shows the motion of components along the 
jet, as wSl as their position angle (PA) with distance 
from the core. We observe that components El-Ela (and 
presumably EO) are ejected with a PA of around —110° 
and they travel towards the south-west direction while 
the other components at 15 GHz are ejected with a PA 
between —117° and —123° and move initially in a less 
southern direction. Gomponents E5 and E 6 present a 
significant change in their velocity vectors at a distance 
of ^3 mas from the core, as observed previously for other 
components and interpreted as the interaction with the 
external medium or a cloud (Gomez et al. 2000, 2008). 
Gomponents dll and dl2-dl2a are ejected with a PA of 
around — 110 °, similar to that observed previously for 
El-Ela. Note that in a precession model components 
El and dl2 would have similar projected velocity vec¬ 
tors while differing in their orientation with respect to 
the observer. Despite the initially different position an¬ 
gles, components move following parallel paths after the 
initial ~ 4-5 mas. Rather than a precession of the whole 
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Fig. 8.— Positions (/eft) and position angles as a function of distance from the core (right) for the model-fit components at 15 (circles) 
and 43 GHz (squares). For clarity errors bars are shown only in the plot of the components’ position. 


jet, we therefore favor the model in which the jet consists 
of a broad funnel through which components - not filling 
the whole jet width - are ejected and travel at differ po¬ 
sition angl es, as supported also by previous observations 
of 3C 120 (Gomez et al. 2011|) and o ther sources in the 
MOJAVE sample (jhister et al.||2U13l). 

Considering that the slower apparent velocities of dll 
and dl 2 are related to the first 7 -ray detections, we 
conclude that these components are most likely ejected 
with a smaller viewing angle. This should increase the 
Doppler factor, leading to the enhanced 7 -ray emission 
measured since the end of 2012 , and a significant increase 
in the total flux of the core and innermost stationary fea¬ 
tures (see Fig. However, we cannot completely rule 
out the possibility that the smaller apparent velocities 
are just due to greater viewing angles, or slowing of the 
components’ velocity. 

It is possible to estimate the required minimum 
Lorentz factor, Tmin = (1 -|- = 6.3, using the 

observed maximum apparent velocity of 6.2c. To min¬ 
imize the required reorientation of the components, we 
can assume that the maximum apparent velocity is ob¬ 
tained for the angle that actually maximizes the apparent 
velocity, given hy 0 = arccos(/3), where /3 is the compo¬ 
nent’s velocity in units of the speed of light. Using our 
previous estimation of Tmin = 6.3, we obtain a viewing 
angle of 0 ~ 9.2° for the maximum apparent velocity of 
6.2c measured for component EO, in agreement with |Ho^ 


vatta et al. (2009). Given that the observed decrease in 
the apparent velocity is probably produced by a decrease 
in the viewing angle, the smallest observed apparent ve¬ 
locity of 4.2c for dll requires 9 ~ 3.6°. These values 
agree with the precession model of Gaproni & Abraham 
(2004), for which the authors estimate T = 6.8 ±0.5 and 
a variation of the viewing angle between (6.3±0.8°) and 
(3.3±0.8°). We should note however that our measured 
change in the apparent velo cities is shifted from the pe- 
riodicity phase predicted by Gaproni & Abraham (20041 
model. 

This change in the orientation from 9 ~ 9.2° to 0 ~ 
3.6° would lead to an increase in the Doppler factor from 
6 ^ 6.2 to (5 ^ 10.9, enhancing the 7 -ray emission above 
the flux detectable by Eermi. 

In summary, all these evidence suggest that the ob¬ 
served 7 -ray emission in 3C 120 depends strongly on the 
orientation of component’s motion with respect to the 


observer, so that only when they are best oriented and 
a new superluminal component pass through the radio 
core a clear 7 -ray detection is obtained. 

Superluminal components are associated with shocks 
moving along the jet, hence the inferred values for the 
velocity - provided there is an estimation of the jet orien¬ 
tation - correspond to the pattern velocity of the shock, 
not the actual flow velocity. The detection of the for¬ 
ward and reserve shocks of the perturbation associated 
with dl 2 and dl 2 a, respectively, allows to obtain a di¬ 
rect calculation of the jet bulk flow velocity. For the 
estimated viewing angle of ~ 3.6° during the ejection of 
components dl 2 and dl 2 a we obtain the corresponding 
pattern velocities of rdi 2 = 6.7 and rdi 2 a= 5 . Hence, we 
can conclude that the jet bulk flow velocity is restricted 
to be 5 < Ej < 6.7. 


6 . 2 . Gamma-ray location and emission mechanism 

Our Hrst 7 -ray detection at the end of 2012 is associ¬ 
ated with the passage of component dll through the mm- 
VLBI core, whose time of ejection is coincident with the 
3 months-bin 7 -ray detection, and is ~ 34 days (0.09 yr) 
after the 15 day s -bin 7 -ray detection. On the other hand, 
Marscher et al. (2002) and Chatterjee et al. (2009) have 
measured a mean delay of 0.18 yr (~ 66 days) between 
X-ray dips and the ejection of new superluminal compo¬ 
nents. Considering that most of the X-ray emission in 
3C 120 originates from the disk-corona system, close to 
the central black hole (BH), we can consider this time 
delay as the distance in time between the BH and the 
mm-VLBI core. A mean apparent motion of ^ 2 mas/yr 
corresponds to a rate of increase in projected separation 
from the core of ~0.24 pc/yr - or slightly smaller if we 
allow for some initial acceleration of the components. 

The measured delay of 34 days between the 15 days-bin 
7 -ray detection and the passage of dll from the mm- 
VLBI core locates the 7 -ray emission upstream of the 
mm-VLBI core, at a distance about half of that between 
the BH and the core. Relativ istic 3D hydrodynamic 
simulations ofjAloy et al. (2003) show that time delays 
stretch signihcantly the observed size of components, so 
that our VLBI observations could detect only some por¬ 
tions of the component, depending on the strength of the 
forward/reverse shock and Doppler factor (i.e., viewing 
angle). Considering also that we have estimated for com¬ 
ponent dll a size of ~0.35 mas (0.2 pc) once it is clearly 
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detached from the core, it is possible that the 15 days- 
bin 7 -ray detection could correspond to the passing of 
the forward section of the dll perturbation through the 
mm-VLBl core, while later on our VLBl images identify 
only the back section , as it is precisely seen in the nu¬ 
merical simulations of Aloy et al. (20031. In this case the 
flaring in y-rays would more closely mark the crossing of 
the dll perturbation through the mm-VLBl core. 

The second y-ray event is instead associated with the 
ejection of components dl2 and dl2a. In this case the 
time of ejection of these two components is also coin¬ 
cident with the 3 months bin y-ray detection, but the 
finer time sampling of the 15-days y-ray light curve con¬ 
strains the y-ray production to ~ 33 days after. Hence 
we conclude that in the second y-ray event the high en¬ 
ergy emission is produced downstream of the mm-VLBl 
core, at a projected distance of ~ 0.13 pc from its posi¬ 
tion - smaller if we consider that the shock responsible 
for the y-ray emission is the slower moving dl2a. 

Reverberation mapping studies of the broad-line- 
region (BLR) in 3C 120 suggest an inclined disk model 
with an extension from 12±7 light days ( ^0.01 pc) to 


28±9 light days (^0.025 pc), from the BH (Grier et al. 


2013 Kollatschny et al. 


2014). These studies also found 


evidence of radial stratihcation in the BLR as well as in¬ 
fall and rotation related to the BH gravity. Therefore the 
BLR extends to about half of the estimated distance be¬ 
tween the BH and the mm-VLBl core, severely limiting 
the external photon field from the BLR at the location 
of the mm-VLBl core. Hence we favor synchrotron self 
Compton as the mechanism for the production of y-ray 
photons in the case of the second y-ray event in late 
2013, although we cannot discard the contribution from 
another external photon field, such the sheath or the ex¬ 
ternal ionized cloud. 

The External Compton process with photons coming 


from the BLR could be instead a possibility for the first 
y-ray event, as our findings point to a y-ray dissipation 
zone between the BH and the radio core, on the limit of 
the extension of the BLR. 

The different orientation of the components within the 
broader jet funnel supports also a model i n which they 
interact wit h the sheath of the jet (see also Gomez et al. 
2000 2008). In this case the observed y-ray activity can 
be produced by Compton s cattering of pho t ons fr om the 


jet sheath, as proposed by Marscher et al. 
p lain the y-ray activit y seen in BKK 1510- 


Grandi et al. 


(2010) to ex- 
( 2012b|) obtain similar results for the FRH 


radio galaxy 3C 111, associating the y-ray activity with 
the ejection of a new radio component from the core, 
and confining the y-ray dissipation region within 0.1 pc 
at a distance of almost 0.3 pc from the black hole. These 
two radio galaxies have other similarities as they are both 
BLRGs and they show a con nection between the radio jet 
and the corona-di sk system (Marscher et al.||2002 Chat- 
terjee et aL]|2011 ). In addition, the apparent motions of 
superluminal cornponents detected in their jets detach 
from th ose of the other ra dio galaxies in the MOJAVE 
sample (Lister et al.||2013). 
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TABLE 3 - VLBA 15 GHz Model-fit Components’ Parameters 


Epoch 

Epoch 

Name 

Flux 

Distance from 

Pos. Angle 

FWHM 

(year) 

(MJD) 


(mjy) 

CO 

(mas) 

(° 



(mas) 

2008.48 

54642.50 

CO 

768 ± 38 







0.31 

± 

0.20 



Cl 

271 ± 14 

0.65 

± 

0.08 

-115.2 

± 

6.7 

0.49 

± 

0.16 



EO 

117 ± 6 

3.51 

± 

0.24 

-112.1 

± 

3.6 

0.79 

± 

0.47 



El 

392 ± 20 

1.50 

± 

0.04 

-109.3 

± 

1.4 

0.32 

± 

0.08 



Ela 

429 ± 21 

1.14 

± 

0.04 

-106.2 

± 

1.7 

0.33 

± 

0.08 

2008.58 

54677.50 

CO 

604 ± 30 







0.22 

± 

0.17 



Cl 

443 ± 22 

0.35 

± 

0.03 

-115.8 

± 

5.3 

0.31 

± 

0.07 



EO 

123 ± 6 

3.77 

± 

0.26 

-111.2 

± 

3.7 

0.88 

± 

0.53 



El 

455 ± 23 

1.75 

± 

0.04 

-109.9 

± 

1.3 

0.38 

± 

0.09 



Ela 

740 ± 37 

1.31 

± 

0.04 

-110.0 

± 

1.7 

0.48 

± 

0.09 

2008.65 

54703.50 

CO 

713 ± 36 







0.22 

± 

0.18 



Cl 

318 ± 16 

0.40 

± 

0.06 

-113.8 

± 

7.9 

0.41 

± 

0.12 



EO 

110 ± 6 

3.92 

± 

0.27 

-110.4 

± 

3.7 

0.87 

± 

0.55 



El 

401 ± 20 

1.81 

± 

0.04 

-109.0 

± 

1.2 

0.35 

± 

0.09 



Ela 

683 ± 34 

1.35 

± 

0.04 

-109.9 

± 

1.6 

0.46 

± 

0.09 

2008.76 

54742.50 

CO 

1084 ± 54 







0.31 

± 

0.12 



Cl 

201 ± 10 

0.46 

± 

0.09 

-116.3 

± 

10.1 

0.46 

± 

0.18 



EO 

127 ± 6 

4.03 

± 

0.24 

-111.9 

± 

3.1 

0.82 

± 

0.47 



El 

586 ± 29 

1.79 

± 

0.05 

-107.8 

± 

1.4 

0.47 

± 

0.10 



Ela 

491 ± 25 

1.25 

± 

0.05 

-110.7 

± 

2.2 

0.45 

± 

0.10 

2008.90 

54796.50 

CO 

492 ± 25 







0.18 

± 

0.09 



Cl 

532 ± 27 

0.45 

± 

0.03 

-122.6 

± 

4.2 

0.33 

± 

0.07 



EO 

122 ± 6 

4.48 

± 

0.22 

-112.4 

± 

2.6 

0.75 

± 

0.43 



El 

223 ± 11 

2.45 

± 

0.06 

-110.3 

± 

1.2 

0.34 

± 

0.12 



S2 

390 ± 20 

1.87 

± 

0.06 

-110.6 

± 

2.6 

0.45 

± 

0.12 



S3 

415 ± 21 

1.25 

± 

0.06 

-112.1 

± 

2.7 

0.49 

± 

0.13 

2009.08 

54861.50 

CO 

395 ± 20 







0.19 

± 

0.10 



Cl 

287 ± 14 

0.49 

± 

0.02 

-123.0 

± 

3.0 

0.14 

± 

0.05 



EO 

134 ± 7 

4.91 

± 

0.24 

-113.7 

± 

2.6 

0.86 

± 

0.48 



El 

118 ± 6 

3.01 

± 

0.11 

-111.6 

± 

1.8 

0.40 

± 

0.21 



E4 

757 ± 38 

0.98 

± 

0.02 

-122.8 

± 

0.9 

0.22 

± 

0.03 



S2 

180 ± 9 

2.22 

± 

0.11 

-111.9 

± 

2.6 

0.52 

± 

0.22 



S3 

418 ± 21 

1.47 

± 

0.04 

-111.7 

± 

3.0 

0.35 

± 

0.08 

2009.23 

54915.50 

CO 

434 ± 22 







0.31 

± 

0.15 



Cl 

176 ± 9 

0.67 

± 

0.02 

-120.6 

± 

1.5 

0.11 

± 

0.05 



EO 

133 ± 7 

5.24 

± 

0.27 

-113.2 

± 

2.7 

0.94 

± 

0.54 



El 

126 ± 6 

3.41 

± 

0.17 

-112.9 

± 

2.7 

0.63 

± 

0.35 



E4 

1385 ± 69 

1.26 

± 

0.02 

-117.7 

± 

0.8 

0.33 

± 

0.05 



S2 

68 ± 3 

2.47 

± 

0.22 

-113.4 

± 

4.8 

0.57 

± 

0.44 



S3 

230 ± 11 

1.81 

± 

0.01 

-108.7 

± 

1.5 

0.11 

± 

0.05 

2009.40 

54979.50 

CO 

393 ± 20 







0.21 

± 

0.09 



Cl 

255 ± 13 

0.35 

± 

0.03 

-117.1 

± 

5.5 

0.23 

± 

0.07 



EO 

165 ± 8 

5.73 

± 

0.32 

-113.6 

± 

3.0 

1.19 

± 

0.63 



El 

139 ± 7 

3.77 

± 

0.28 

-110.9 

± 

3.9 

0.98 

± 

0.56 



E4 

1489 ± 74 

1.61 

± 

0.02 

-116.0 

± 

0.6 

0.32 

± 

0.04 



E4a 

698 ± 35 

1.22 

± 

0.03 

-115.9 

± 

1.4 

0.37 

± 

0.06 



S3 

106 ± 5 

2.25 

± 

0.12 

-108.2 

± 

2.8 

0.43 

± 

0.24 

2009.51 

55017.50 

CO 

285 ± 14 







0.24 

± 

0.09 



Cl 

287 ± 14 

0.48 

± 

0.05 

-116.4 

± 

5.7 

0.34 

± 

0.10 



EO 

140 ± 7 

6.02 

± 

0.35 

-113.8 

± 

3.1 

1.19 

± 

0.70 



El 

136 ± 7 

4.08 

± 

0.29 

-111.1 

± 

3.8 

1.01 

± 

0.58 



E4 

756 ± 38 

1.80 

± 

0.03 

-114.1 

± 

0.8 

0.33 

± 

0.05 



E4a 

896 ± 45 

1.48 

± 

0.03 

-114.6 

± 

1.1 

0.38 

± 

0.06 



S3 

107 ± 5 

2.26 

± 

0.23 

-114.5 

± 

5.4 

0.73 

± 

0.45 

2009.63 

55062.50 

CO 

470 ± 24 







0.26 

± 

0.08 



EO 

133 ± 7 

6.27 

± 

0.40 

-113.9 

± 

3.4 

1.30 

± 

0.79 



El 

137 ± 7 

4.29 

± 

0.33 

-110.7 

± 

4.5 

1.13 

± 

0.67 



E4 

839 ± 42 

1.93 

± 

0.03 

-112.6 

± 

0.9 

0.41 

± 

0.07 
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E4a 

E5 

635 ± 32 
266 ± 13 

1.55 ± 0.03 
0.70 ± 0.08 

-114.6 ± 1.2 
-116.7 ± 3.1 

0.36 ± 0.07 
0.49 ± 0.16 

2009.94 

55175.50 

CO 

463 ± 23 



0.20 ± 0.08 



Cl 

179 ± 9 

0.46 ± 0.08 

-119.5 ± 9.8 

0.40 ± 0.16 



EO 

78 ± 4 

7.13 ± 0.53 

-114.8 ± 4.0 

1.26 ± 1.06 



El 

117 ± 6 

5.30 ± 0.40 

-111.3 ± 4.0 

1.22 ± 0.79 



E4 

849 ± 42 

2.31 ± 0.07 

-114.1 ± 1.5 

0.74 ± 0.13 



E5 

523 ± 26 

1.45 ± 0.06 

-113.3 ± 2.3 

0.56 ± 0.13 

2010.11 

55238.50 

CO 

353 ± 18 



0.16 ± 0.10 



Cl 

204 ± 10 

0.46 ± 0.03 

-118.5 ± 3.9 

0.20 ± 0.07 



EO 

130 ± 7 

7.61 ± 0.60 

-115.5 ± 4.4 

1.86 ± 1.20 



El 

109 ± 5 

5.44 ± 0.47 

-111.3 ± 3.9 

1.34 ± 0.94 



E4 

602 ± 30 

2.79 ± 0.11 

-114.4 ± 2.1 

0.93 ± 0.22 



E5 

628 ± 31 

1.78 ± 0.05 

-113.3 ± 1.3 

0.46 ± 0.09 



S7 

151 ± 8 

1.11 ± 0.13 

-114.8 ± 6.3 

0.54 ± 0.26 

2010.19 

55265.50 

CO 

444 ± 22 



0.24 ± 0.08 



Cl 

136 ± 7 

0.53 ± 0.06 

-118.2 ± 5.8 

0.25 ± 0.11 



EO 

87 ± 4 

7.80 ± 0.50 

-116.5 ± 3.5 

1.26 ± 0.99 



El 

95 ± 5 

6.00 ± 0.43 

-111.2 ± 4.5 

1.17 ± 0.86 



E4 

392 ± 20 

2.97 ± 0.13 

-113.1 ± 2.4 

0.90 ± 0.27 



E5 

540 ± 27 

1.91 ± 0.05 

-112.5 ± 1.3 

0.44 ± 0.09 



S7 

279 ± 14 

1.22 ± 0.05 

-115.9 ± 2.3 

0.34 ± 0.10 

2010.53 

55389.50 

CO 

519 ± 26 



0.24 ± 0.06 



Cl 

489 ± 24 

0.36 ± 0.03 

-130.6 ± 4.2 

0.26 ± 0.05 



EO 

75 ± 4 

8.83 ± 0.60 

-116.4 ± 4.4 

1.58 ± 1.20 



El 

77 ± 4 

6.89 ± 0.37 

-111.8 ± 3.8 

0.92 ± 0.74 



E4 

204 ± 10 

3.93 ± 0.18 

-114.1 ± 2.4 

0.81 ± 0.35 



E5 

197 ± 10 

2.65 ± 0.17 

-111.1 ± 3.3 

0.77 ± 0.34 



E6 

284 ± 14 

1.81 ± 0.06 

-113.5 ± 1.8 

0.39 ± 0.12 



S7 

180 ± 9 

1.31 ± 0.08 

-117.4 ± 3.2 

0.38 ± 0.15 

2010.65 

55435.50 

CO 

345 ± 17 



0.21 ± 0.07 



Cl 

401 ± 20 

0.48 ± 0.02 

-124.9 ± 2.9 

0.22 ± 0.05 



EO 

45 ± 2 

9.83 ± 0.29 

-117.8 ± 1.6 

0.58 ± 0.59 



El 

121 ± 6 

7.38 ± 0.60 

-113.4 ± 2.8 

1.76 ± 1.20 



E4 

131 ± 7 

4.38 ± 0.18 

-115.5 ± 2.3 

0.67 ± 0.37 



E5 

174 ± 9 

3.28 ± 0.26 

-109.3 ± 4.1 

1.04 ± 0.52 



E6 

297 ± 15 

2.04 ± 0.09 

-114.5 ± 2.4 

0.56 ± 0.18 



S7 

108 ± 5 

1.29 ± 0.16 

-116.2 ± 6.5 

0.53 ± 0.31 

2010.74 

55468.50 

CO 

439 ± 22 



0.23 ± 0.08 



Cl 

297 ± 15 

0.60 ± 0.05 

-120.3 ± 4.3 

0.32 ± 0.09 



EO 

62 ± 3 

10.10 ± 0.25 

-118.8 ± 1.4 

0.59 ± 0.50 



El 

104 ± 5 

7.46 ± 0.41 

-112.7 ± 4.4 

1.17 ± 0.82 



E4 

122 ± 6 

4.54 ± 0.18 

-116.3 ± 2.1 

0.64 ± 0.36 



E5 

160 ± 8 

3.62 ± 0.20 

-108.0 ± 2.8 

0.79 ± 0.39 



E6 

235 ± 12 

2.21 ± 0.12 

-112.3 ± 2.9 

0.64 ± 0.24 



S7 

131 ± 7 

1.34 ± 0.11 

-112.6 ± 4.4 

0.45 ± 0.23 

2010.84 

55504.50 

CO 

586 ± 29 



0.24 ± 0.07 



EO 

63 ± 3 

10.39 ± 0.17 

-118.9 ± 0.9 

0.44 ± 0.34 



El 

103 ± 5 

8.07 ± 0.60 

-114.2 ± 2.9 

1.93 ± 1.20 



E4 

152 ± 8 

4.74 ± 0.20 

-115.2 ± 2.3 

0.78 ± 0.40 



E5 

116 ± 6 

3.85 ± 0.16 

-107.7 ± 2.1 

0.56 ± 0.32 



E6 

208 ± 10 

2.41 ± 0.15 

-111.9 ± 3.3 

0.72 ± 0.30 



E8 

231 ± 12 

0.75 ± 0.10 

-120.7 ± 4.1 

0.53 ± 0.19 



S7 

170 ± 9 

1.48 ± 0.09 

-116.2 ± 3.5 

0.44 ± 0.19 

2010.91 

55529.50 

CO 

529 ± 26 



0.20 ± 0.06 



Cl 

104 ± 5 

0.39 ± 0.05 

-128.8 ± 7.9 

0.21 ± 0.11 



EO 

62 ± 3 

10.60 ± 0.23 

-118.7 ± 1.2 

0.57 ± 0.47 



El 

105 ± 5 

8.23 ± 0.60 

-113.7 ± 5.0 

1.93 ± 1.20 



E4 

138 ± 7 

4.88 ± 0.24 

-115.8 ± 2.6 

0.85 ± 0.47 



E5 

97 ± 5 

4.01 ± 0.18 

-107.5 ± 2.3 

0.57 ± 0.36 



E6 

139 ± 7 

2.62 ± 0.21 

-111.6 ± 4.2 

0.77 ± 0.41 



E8 

168 ± 8 

1.07 ± 0.09 

-118.2 ± 4.4 

0.41 ± 0.18 



S7 

154 ± 8 

1.61 ± 0.08 

-115.2 ± 3.4 

0.37 ± 0.16 



14 


2010.98 

55554.50 

CO 

563 ± 28 



0.17 ± 0.05 



Cl 

231 ± 12 

0.34 ± 0.03 

-123.3 ± 4.3 

0.18 ± 0.05 



EO 

38 ± 2 

10.64 ± 0.18 

-116.9 ± 1.0 

0.35 ± 0.36 



El 

102 ± 5 

8.41 ± 0.48 

-114.5 ± 4.8 

1.33 ± 0.96 



E4 

132 ± 7 

5.11 ± 0.20 

-116.9 ± 2.1 

0.71 ± 0.39 



E5 

132 ± 7 

4.24 ± 0.21 

-108.6 ± 2.5 

0.76 ± 0.42 



E6 

123 ± 6 

2.84 ± 0.24 

-112.7 ± 4.4 

0.81 ± 0.48 



E8 

331 ± 17 

1.29 ± 0.06 

-117.7 ± 2.7 

0.44 ± 0.13 



S7 

100 ± 5 

1.82 ± 0.05 

-116.3 ± 3.2 

0.19 ± 0.09 

2011.14 

55612.50 

CO 

814 ± 41 



0.18 ± 0.05 



Cl 

187 ± 9 

0.42 ± 0.05 

-117.9 ± 7.2 

0.29 ± 0.11 



EO 

65 ± 3 

11.21 ± 0.23 

-118.9 ± 1.1 

0.56 ± 0.45 



El 

102 ± 5 

8.94 ± 0.60 

-114.7 ± 3.1 

1.90 ± 1.20 



E4 

145 ± 7 

5.38 ± 0.27 

-114.5 ± 2.7 

0.98 ± 0.54 



E5 

53 ± 3 

4.50 ± 0.25 

-107.5 ± 2.9 

0.56 ± 0.51 



E6 

101 ± 5 

3.36 ± 0.30 

-107.6 ± 4.6 

0.90 ± 0.60 



E8 

410 ± 20 

1.54 ± 0.06 

-115.2 ± 2.1 

0.48 ± 0.12 



S7 

93 ± 5 

2.08 ± 0.13 

-115.6 ± 3.2 

0.41 ± 0.25 

2011.16 

55619.50 

CO 

1017 ± 51 



0.18 ± 0.05 



Cl 

225 ± 11 

0.40 ± 0.05 

-123.2 ± 6.9 

0.30 ± 0.10 



EO 

76 ± 4 

11.27 ± 0.21 

-118.9 ± 1.0 

0.59 ± 0.43 



El 

121 ± 6 

8.93 ± 0.60 

-114.7 ± 4.4 

1.85 ± 1.20 



E4 

178 ± 9 

5.41 ± 0.25 

-113.8 ± 2.5 

1.01 ± 0.50 



E5 

52 ± 3 

4.52 ± 0.26 

-108.2 ± 2.9 

0.57 ± 0.52 



E6 

112 ± 6 

3.43 ± 0.24 

-108.3 ± 3.5 

0.77 ± 0.47 



E8 

463 ± 23 

1.57 ± 0.06 

-115.3 ± 2.0 

0.49 ± 0.12 

2011.17 

55625.50 

CO 

954 ± 48 



0.16 ± 0.06 



Cl 

386 ± 19 

0.32 ± 0.02 

-124.8 ± 4.0 

0.20 ± 0.05 



EO 

80 ± 4 

11.29 ± 0.21 

-118.4 ± 1.0 

0.59 ± 0.42 



El 

107 ± 5 

9.02 ± 0.60 

-114.6 ± 3.9 

1.97 ± 1.20 



E4 

186 ± 9 

5.45 ± 0.28 

-113.8 ± 2.7 

1.13 ± 0.55 



E5 

44 ± 2 

4.68 ± 0.14 

-108.8 ± 1.5 

0.30 ± 0.27 



E6 

162 ± 8 

3.35 ± 0.32 

-108.9 ± 4.9 

1.19 ± 0.64 



E8 

385 ± 19 

1.78 ± 0.05 

-114.3 ± 1.5 

0.39 ± 0.10 



S8 

199 ± 10 

1.31 ± 0.08 

-118.1 ± 3.2 

0.40 ± 0.15 

2011.38 

55702.50 

CO 

689 ± 34 



0.17 ± 0.09 



Cl 

343 ± 17 

0.49 ± 0.05 

-116.7 ± 5.1 

0.34 ± 0.09 



EO 

64 ± 3 

11.70 ± 0.29 

-117.9 ± 1.4 

0.69 ± 0.58 



El 

61 ± 3 

9.51 ± 0.60 

-116.5 ± 4.5 

1.37 ± 1.20 



E4 

85 ± 5 

6.02 ± 0.35 

-112.4 ± 3.1 

1.06 ± 0.70 



E5 

0.08 ± 4 

5.09 ± 0.42 

-114.2 ± 4.4 

1.08 ± 0.84 



E6 

110 ± 6 

3.75 ± 0.20 

-106.3 ± 2.7 

0.67 ± 0.41 



E8 

301 ± 15 

2.08 ± 0.06 

-112.8 ± 1.5 

0.39 ± 0.12 



S8 

93 ± 5 

1.41 ± 0.10 

-120.4 ± 3.9 

0.34 ± 0.20 

2011.56 

55766.50 

CO 

483 ± 24 



0.24 ± 0.10 



EO 

59 ± 3 

11.96 ± 0.34 

-117.6 ± 1.5 

0.75 ± 0.67 



El 

68 ± 3 

9.94 ± 0.60 

-115.5 ± 3.9 

1.93 ± 1.20 



E4 

86 ± 4 

6.54 ± 0.38 

-112.3 ± 3.1 

1.01 ± 0.77 



E5 

67 ± 3 

5.50 ± 0.39 

-113.6 ± 3.8 

0.91 ± 0.79 



E6 

144 ± 7 

4.05 ± 0.26 

-109.1 ± 3.3 

0.93 ± 0.51 



E8 

162 ± 8 

2.40 ± 0.08 

-112.3 ± 1.7 

0.36 ± 0.15 



E9 

214 ± 11 

0.71 ± 0.09 

-117.8 ± 4.0 

0.49 ± 0.19 



S8 

88 ± 4 

1.71 ± 0.18 

-115.9 ± 5.6 

0.54 ± 0.35 

2011.65 

55799.50 

CO 

505 ± 25 



0.12 ± 0.08 



Cl 

136 ± 7 

0.45 ± 0.07 

-121.9 ± 9.1 

0.32 ± 0.15 



EO 

55 ± 3 

12.26 ± 0.39 

-117.7 ± 1.7 

0.82 ± 0.78 



El 

68 ± 3 

10.38 ± 0.60 

-115.4 ± 5.2 

1.95 ± 1.20 



E4 

75 ± 4 

6.92 ± 0.42 

-112.3 ± 3.2 

1.01 ± 0.83 



E5 

73 ± 4 

5.90 ± 0.41 

-113.8 ± 3.7 

0.98 ± 0.82 



E6 

138 ± 7 

4.34 ± 0.25 

-110.8 ± 3.0 

0.89 ± 0.49 



E8 

190 ± 10 

2.55 ± 0.15 

-112.9 ± 3.2 

0.70 ± 0.31 



E9 

161 ± 8 

1.18 ± 0.13 

-116.5 ± 4.0 

0.57 ± 0.27 

2011.75 

55837.50 

CO 

625 ± 31 



0.14 ± 0.08 



Cl 

146 ± 7 

0.37 ± 0.06 

-129.8 ± 9.1 

0.27 ± 0.12 



15 




EO 

47 ± 2 

12.49 ± 0.45 

-117.8 ± 1.9 

0.86 ± 0.91 



El 

72 ± 4 

10.62 ± 0.60 

-115.4 ± 5.1 

2.19 ± 1.20 



E4 

69 ± 3 

7.14 ± 0.46 

-111.5 ± 3.4 

1.06 ± 0.92 



E5 

63 ± 3 

6.13 ± 0.36 

-115.0 ± 3.2 

0.82 ± 0.72 



E6 

158 ± 8 

4.57 ± 0.26 

-111.0 ± 3.0 

0.99 ± 0.52 



E8 

108 ± 5 

2.82 ± 0.11 

-111.3 ± 2.0 

0.39 ± 0.21 



E9 

195 ± 10 

1.27 ± 0.08 

-117.8 ± 3.3 

0.39 ± 0.15 

2011.94 

55907.50 

CO 

420 ± 21 



0.17 ± 0.08 



Cl 

143 ± 7 

0.45 ± 0.06 

-124.8 ± 7.9 

0.29 ± 0.13 



EO 

38 ± 2 

12.78 ± 0.60 

-117.8 ± 3.3 

1.20 ± 1.20 



El 

55 ± 3 

10.71 ± 0.60 

-114.7 ± 5.5 

2.28 ± 1.20 



E4 

82 ± 4 

7.11 ± 0.49 

-113.3 ± 3.7 

1.21 ± 0.98 



E6 

140 ± 7 

5.13 ± 0.30 

-112.1 ± 3.1 

1.04 ± 0.60 



E8 

88 ± 4 

3.21 ± 0.20 

-110.6 ± 3.2 

0.59 ± 0.39 



E9 

256 ± 13 

1.65 ± 0.10 

-116.3 ± 3.3 

0.57 ± 0.20 

2011.99 

55924.50 

CO 

444 ± 22 



0.16 ± 0.07 



Cl 

170 ± 8 

0.35 ± 0.06 

-109.7 ± 9.2 

0.31 ± 0.13 



EO 

36 ± 2 

13.03 ± 0.60 

-117.4 ± 3.4 

1.22 ± 1.20 



El 

55 ± 3 

11.15 ± 0.60 

-115.5 ± 6.7 

1.72 ± 1.20 



E4 

88 ± 4 

7.17 ± 0.53 

-112.9 ± 3.9 

1.33 ± 1.06 



E6 

146 ± 7 

5.23 ± 0.30 

-111.7 ± 7.1 

1.08 ± 0.61 



E8 

101 ± 5 

3.32 ± 0.20 

-110.1 ± 3.2 

0.64 ± 0.40 



E9 

154 ± 13 

1.71 ± 0.11 

-115.2 ± 3.5 

0.61 ± 0.22 

2012.04 

55940.50 

CO 

494 ± 25 



0.16 ± 0.07 



Cl 

144 ± 7 

0.48 ± 0.08 

-121.6 ± 9.3 

0.35 ± 0.16 



EO 

38 ± 1 

12.95 ± 0.60 

-117.4 ± 3.8 

1.37 ± 1.20 



El 

49 ± 2 

11.17 ± 0.60 

-115.1 ± 4.6 

1.94 ± 1.20 



E4 

87 ± 4 

7.31 ± 0.60 

-113.1 ± 4.4 

1.49 ± 1.20 



E6 

135 ± 7 

5.39 ± 0.28 

-111.9 ± 3.1 

0.97 ± 0.56 



E8 

116 ± 6 

3.44 ± 0.19 

-111.1 ± 2.9 

0.66 ± 0.39 



E9 

228 ± 11 

1.80 ± 0.13 

-115.6 ± 4.0 

0.68 ± 0.26 

2012.17 

55990.50 

CO 

559 ± 28 



0.16 ± 0.06 



Cl 

111 ± 6 

0.70 ± 0.14 

-121.6 ± 28.6 

0.50 ± 0.29 



EO 

24 ± 2 

12.97 ± 0.41 

-118.3 ± 1.7 

0.56 ± 0.82 



El 

48 ± 2 

11.89 ± 0.60 

-114.6 ± 5.7 

2.25 ± 1.20 



E4 

70 ± 4 

7.75 ± 0.60 

-112.8 ± 4.9 

1.53 ± 1.20 



E6 

128 ± 6 

5.77 ± 0.31 

-111.9 ± 2.7 

1.03 ± 0.62 



E8 

199 ± 10 

3.69 ± 0.11 

-110.2 ± 1.5 

0.53 ± 0.21 



E9 

145 ± 7 

2.02 ± 0.19 

-116.1 ± 5.2 

0.74 ± 0.39 

2012.33 

56046.50 

CO 

560 ± 28 



0.17 ± 0.06 



Cl 

140 ± 7 

0.37 ± 0.06 

-123.9 ± 26.1 

0.29 ± 0.13 



El 

75 ± 4 

12.84 ± 0.60 

-116.1 ± 6.4 

1.93 ± 1.20 



E4 

66 ± 3 

8.11 ± 0.60 

-113.4 ± 5.0 

1.58 ± 1.20 



E6 

123 ± 6 

6.25 ± 0.33 

-112.4 ± 3.8 

1.06 ± 0.66 



E8 

279 ± 14 

3.98 ± 0.09 

-110.8 ± 1.2 

0.54 ± 0.18 



E9 

94 ± 5 

2.47 ± 0.27 

-115.1 ± 5.9 

0.79 ± 0.54 



S9 

68 ± 3 

1.36 ± 0.28 

-120.0 ± 8.0 

0.68 ± 0.55 

2012.39 

56071.50 

CO 

509 ± 25 



0.13 ± 0.05 



Cl 

162 ± 8 

0.32 ± 0.05 

-125.6 ± 8.0 

0.24 ± 0.09 



El 

63 ± 3 

12.96 ± 0.60 

-116.4 ± 3.8 

1.67 ± 1.20 



E4 

42 ± 2 

8.62 ± 0.60 

-113.4 ± 5.1 

1.34 ± 1.20 



E6 

138 ± 7 

6.42 ± 0.53 

-112.4 ± 2.8 

1.67 ± 1.05 



E8 

289 ± 14 

4.14 ± 0.08 

-111.1 ± 1.0 

0.51 ± 0.16 



E9 

80 ± 4 

2.71 ± 0.21 

-115.9 ± 4.2 

0.59 ± 0.42 



S9 

73 ± 4 

1.44 ± 0.30 

-118.3 ± 8.0 

0.75 ± 0.59 

2012.53 

56120.50 

CO 

548 ± 27 



0.15 ± 0.06 



Cl 

166 ± 8 

0.27 ± 0.02 

-118.0 ± 4.6 

0.13 ± 0.05 



C2 

50 ± 3 

0.65 ± 0.13 

-120.2 ± 6.0 

0.31 ± 0.26 



El 

57 ± 3 

13.31 ± 0.60 

-116.3 ± 3.5 

1.78 ± 1.20 



E4 

47 ± 2 

8.68 ± 0.60 

-113.5 ± 5.1 

1.40 ± 1.20 



E6 

99 ± 5 

6.78 ± 0.42 

-112.4 ± 4.3 

1.17 ± 0.84 



E8 

263 ± 13 

4.46 ± 0.13 

-111.8 ± 1.5 

0.70 ± 0.25 



E9 

60 ± 3 

2.99 ± 0.50 

-113.1 ± 7.0 

1.05 ± 0.10 



S9 

45 ± 2 

1.68 ± 0.52 

-118.7 ± 9.5 

0.95 ± 1.04 
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2012.59 

56142.50 

CO 

Cl 

C2 

El 

E4 

E6 

E8 

E9 

S9 

490 ± 25 
319 ± 16 

49 ± 2 

50 ± 2 

43 ± 2 

93 ± 5 

238 ± 12 

69 ± 3 

46 ± 2 

0.23 ± 0.01 
0.65 ± 0.04 

13.55 ± 0.60 
8.90 ± 0.60 
6.99 ± 0.42 
4.66 ± 0.13 
3.34 ± 0.57 

1.55 ± 0.56 

-115.0 ± 2.5 
-116.9 ± 3.1 

-115.7 ± 3.9 
-113.5 ± 5.9 
-112.5 ± 3.3 
-111.9 ± 1.5 
-112.5 ± 7.0 
-118.7 ± 9.5 

< 0.05 

0.09 ± 0.02 
0.11 ± 0.07 
1.65 ± 1.20 
1.56 ± 1.20 
1.13 ± 0.83 
0.70 ± 0.27 
1.26 ± 1.15 
1.02 ± 1.13 

2012.67 

56172.50 

CO 

582 ± 29 



< 0.05 



Cl 

424 ± 21 

0.22 ± 0.01 

-117.5 ± 1.6 

0.08 ± 0.01 



C2 

63 ± 3 

0.64 ± 0.06 

-115.9 ± 5.5 

0.20 ± 0.13 



El 

65 ± 3 

13.69 ± 0.60 

-115.5 ± 3.8 

2.29 ± 1.20 



E4 

40 ± 2 

9.35 ± 0.60 

-113.8 ± 6.0 

1.60 ± 1.20 



E6 

100 ± 5 

7.24 ± 0.50 

-112.5 ± 3.2 

1.35 ± 0.99 



E8 

208 ± 10 

4.91 ± 0.17 

-112.1 ± 1.8 

0.79 ± 0.34 



E9 

68 ± 3 

3.59 ± 0.58 

-112.4 ± 7.0 

1.27 ± 1.16 



S9 

54 ± 3 

1.53 ± 0.40 

-120.3 ± 9.6 

0.83 ± 0.81 

2012.84 

56233.50 

CO 

641 ± 32 



0.10 ± 0.07 



Cl 

478 ± 24 

0.23 ± 0.01 

-113.1 ± 2.2 

0.11 ± 0.019 



C2 

70 ± 3 

0.59 ± 0.05 

-113.2 ± 4.1 

0.16 ± 0.09 



El 

57 ± 3 

14.17 ± 0.60 

-115.0 ± 4.8 

2.30 ± 1.20 



E4 

34 ± 2 

9.76 ± 0.60 

-113.6 ± 6.2 

1.56 ± 1.20 



E6 

74 ± 4 

7.70 ± 0.60 

-112.2 ± 3.6 

1.40 ± 1.20 



E8 

181 ± 9 

5.32 ± 0.22 

-112.5 ± 2.2 

0.92 ± 0.44 



E9 

56 ± 3 

3.87 ± 0.57 

-112.9 ± 7.0 

1.14 ± 1.15 



S9 

43 ± 2 

1.44 ± 0.48 

-119.6 ± 9.7 

0.87 ± 0.97 

2012.91 

56259.50 

CO 

491 ± 25 



< 0.05 



Cl 

590 ± 29 

0.23 ± 0.01 

-116.5 ± 0.6 

0.05 ± 0.05 



C2 

123 ± 6 

0.55 ± 0.03 

-114.8 ± 3.2 

0.16 ± 0.07 



El 

40 ± 2 

14.29 ± 0.60 

-115.7 ± 4.9 

1.57 ± 1.20 



E4 

26 ± 1 

9.97 ± 0.60 

-113.6 ± 4.6 

1.09 ± 1.20 



E6 

64 ± 3 

8.01 ± 0.60 

-112.3 ± 4.3 

1.42 ± 1.20 



E8 

178 ± 9 

5.53 ± 0.24 

-112.8 ± 2.3 

0.99 ± 0.49 



E9 

52 ± 3 

4.01 ± 0.54 

-113.9 ± 5.7 

1.04 ± 1.07 



S9 

43 ± 2 

1.31 ± 0.43 

-116.4 ± 9.0 

0.79 ± 0.87 

2012.98 

56284.50 

CO 

633 ± 32 



0.06 ± 0.05 



Cl 

535 ± 27 

0.26 ± 0.01 

-118.0 ± 2.0 

< 0.05 



C2 

82 ± 4 

0.70 ± 0.01 

-116.7 ± 3.2 

< 0.05 



El 

45 ± 2 

14.51 ± 0.60 

-114.9 ± 3.9 

1.85 ± 1.20 



E4 

33 ± 2 

10.52 ± 0.60 

-113.4 ± 5.8 

1.56 ± 1.20 



E6 

66 ± 3 

8.01 ± 0.60 

-112.7 ± 4.6 

1.39 ± 1.20 



E8 

176 ± 9 

5.63 ± 0.25 

-112.8 ± 2.3 

1.01 ± 0.49 



E9 

61 ± 3 

4.11 ± 0.44 

-113.7 ± 5.9 

0.95 ± 0.88 



S9 

27 ± 1 

1.59 ± 0.24 

-115.1 ± 8.1 

0.38 ± 0.48 

2013.05 

56313.50 

CO 

665 ± 33 



< 0.05 



Cl 

575 ± 29 

0.23 ± 0.01 

-117.6 ± 2.5 

0.08 ± 0.05 



C2 

109 ± 5 

0.60 ± 0.07 

-115.5 ± 5.8 

0.26 ± 0.13 



El 

55 ± 3 

14.64 ± 0.60 

-114.9 ± 4.4 

2.85 ± 1.20 



E4 

30 ± 2 

10.43 ± 0.60 

-113.9 ± 7.1 

1.73 ± 1.20 



E6 

60 ± 3 

8.38 ± 0.60 

-112.09 ± 4.4 

1.43 ± 1.20 



E8 

170 ± 9 

5.92 ± 0.26 

-113.2 ± 2.3 

1.03 ± 0.52 



E9 

70 ± 4 

4.41 ± 0.50 

-112.2 ± 5.7 

1.13 ± 0.99 



S9 

46 ± 2 

1.40 ± 0.44 

-117.9 ± 9.1 

0.83 ± 0.89 

2013.11 

56333.50 

CO 

822 ± 41 



0.11 ± 0.07 



Cl 

594 ± 30 

0.24 ± 0.01 

-117.6 ± 2.3 

< 0.05 



C2 

73 ± 4 

0.56 ± 0.01 

-114.9 ± 1.3 

0.06 ± 0.05 



El 

41 ± 2 

14.27 ± 0.60 

-114.5 ± 6.4 

2.67 ± 1.20 



E6 

77 ± 4 

8.61 ± 0.60 

-111.8 ± 4.6 

1.61 ± 1.20 



E8 

78 ± 4 

6.31 ± 0.20 

-112.4 ± 1.7 

0.57 ± 0.41 



E8a 

62 ± 3 

5.64 ± 0.20 

-114.6 ± 1.9 

0.49 ± 0.40 



E9 

87 ± 4 

4.52 ± 0.49 

-111.6 ± 5.1 

1.25 ± 0.98 



S9 

41 ± 2 

1.07 ± 0.21 

-117.3 ± 9.2 

0.42 ± 0.42 
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2013.16 

56351.50 

CO 

Cl 

C2 

El 

E6 

E8 

E8a 

E9 

S9 

922 ± 46 
642 ± 32 

88 ± 4 

40 ± 2 

71 ± 4 

84 ± 4 

62 ± 3 

82 ± 4 

53 ± 3 

0.24 ± 0.01 

0.55 ± 0.01 
14.47 ± 0.60 

8.85 ± 0.60 
6.44 ± 0.21 
5.72 ± 0.20 
4.65 ± 0.45 
1.12 ± 0.28 

-121.6 ± 2.1 
-116.5 ± 1.3 
-114.2 ± 7.2 
-112.6 ± 4.4 
-112.5 ± 1.7 
-115.1 ± 1.9 
-111.9 ± 6.1 
-115.4 ± 9.5 

0.10 ± 0.05 

< 0.05 

< 0.05 

2.60 ± 1.20 
1.86 ± 1.20 
0.61 ± 0.42 
0.49 ± 0.39 
1.13 ± 0.90 
0.60 ± 0.55 

2013.24 

56382.50 

CO 

718 ± 36 



< 0.05 



Cl 

1123 ± 56 

0.25 ± 0.01 

-116.9 ± 2.0 

0.09 ± 0.05 



C2 

114 ± 6 

0.60 ± 0.01 

-119.4 ± 2.9 

0.07 ± 0.02 



El 

28 ± 1 

14.67 ± 0.60 

-113.4 ± 7.0 

3.01 ± 1.20 



E6 

53 ± 3 

9.20 ± 0.60 

-111.9 ± 5.4 

1.75 ± 1.20 



E8 

75 ± 4 

6.76 ± 0.22 

-112.3 ± 1.7 

0.59 ± 0.44 



E8a 

70 ± 4 

5.88 ± 0.20 

-114.7 ± 1.8 

0.53 ± 0.39 



E9 

76 ± 4 

4.76 ± 0.55 

-112.2 ± 5.3 

1.28 ± 1.10 



S9 

44 ± 2 

1.43 ± 0.31 

-117.4 ± 9.5 

0.61 ± 0.63 

2013.42 

56445.50 

CO 

530 ± 27 



0.13 ± 0.06 



Cl 

626 ± 31 

0.27 ± 0.01 

-114.9 ± 2.2 

0.14 ± 0.05 



C2 

369 ± 18 

0.63 ± 0.02 

-113.4 ± 1.5 

0.15 ± 0.03 



El 

29 ± 1 

15.27 ± 0.60 

-114.5 ± 7.2 

2.38 ± 1.20 



E6 

55 ± 3 

9.92 ± 0.60 

-112.7 ± 5.7 

2.18 ± 1.20 



E8 

67 ± 3 

7.16 ± 0.40 

-112.2 ± 2.9 

0.92 ± 0.79 



E8a 

68 ± 3 

6.24 ± 0.15 

-114.4 ± 1.3 

0.40 ± 0.30 



E9 

80 ± 4 

5.17 ± 0.51 

-113.3 ± 5.4 

1.24 ± 1.02 



S9 

64 ± 3 

1.37 ± 0.37 

-120.1 ± 9.4 

0.84 ± 0.74 

2013.51 

56481.50 

CO 

687 ± 34 



0.16 ± 0.06 



Cl 

457 ± 23 

0.28 ± 0.02 

-113.5 ± 3.3 

0.18 ± 0.05 



C2 

235 ± 12 

0.79 ± 0.04 

-114.9 ± 2.9 

0.27 ± 0.09 



El 

17 ± 1 

15.56 ± 0.60 

-116.0 ± 7.6 

1.94 ± 1.20 



E6 

26 ± 1 

10.14 ± 0.60 

-112.2 ± 6.7 

1.89 ± 1.20 



E8 

46 ± 3 

7.46 ± 0.48 

-112.3 ± 3.4 

0.89 ± 0.96 



E8a 

70 ± 4 

6.44 ± 0.18 

-113.9 ± 1.5 

0.48 ± 0.35 



E9 

60 ± 3 

5.31 ± 0.55 

-113.5 ± 5.9 

1.14 ± 1.11 



S9 

40 ± 2 

1.50 ± 0.45 

-119.5 ± 9.4 

0.79 ± 0.91 

2013.57 

56503.50 

CO 

694 ± 35 



0.10 ± 0.07 



Cl 

525 ± 26 

0.30 ± 0.01 

-113.6 ± 2.0 

0.13 ± 0.02 



C2 

139 ± 7 

0.72 ± 0.03 

-114.9 ± 2.4 

0.16 ± 0.06 



El 

19 ± 1 

15.56 ± 0.60 

-115.7 ± 4.5 

1.32 ± 1.20 



E6 

31 ± 2 

10.67 ± 0.60 

-112.0 ± 8.7 

1.75 ± 1.20 



E8 

47 ± 2 

7.70 ± 0.49 

-112.8 ± 3.3 

0.91 ± 0.97 



E8a 

76 ± 4 

6.62 ± 0.18 

-113.7 ± 1.4 

0.49 ± 0.35 



E9 

70 ± 4 

5.44 ± 0.60 

-113.9 ± 5.0 

1.32 ± 1.19 



S9 

115 ± 6 

1.07 ± 0.06 

-116.1 ± 7.0 

0.26 ± 0.13 

2013.63 

56524.50 

CO 

752 ± 38 



0.18 ± 0.07 



Cl 

340 ± 17 

0.35 ± 0.03 

-111.4 ± 4.0 

0.22 ± 0.05 



C2 

158 ± 8 

0.98 ± 0.08 

-118.3 ± 4.6 

0.38 ± 0.16 



El 

11 ± 1 

16.46 ± 0.58 

-119.2 ± 1.9 

0.51 ± 1.15 



E6 

34 ± 2 

11.63 ± 1.88 

-113.4 ± 6.8 

2.39 ± 3.76 



E8 

31 ± 2 

7.89 ± 0.43 

-110.4 ± 3.0 

0.66 ± 0.86 



E8a 

80 ± 4 

6.72 ± 0.21 

-112.7 ± 1.7 

0.59 ± 0.42 



E9 

54 ± 3 

5.50 ± 0.51 

-114.1 ± 5.3 

1.02 ± 1.02 



S9 

19 ± 1 

1.71 ± 0.51 

-117.8 ± 9.0 

0.60 ± 1.01 


TABLE 4 - VLBA 43 GHz Model-fit Components’ Parameters 


Epoch 

(year) 

Epoch 

(MJD) 

Name 

Flux 

(mjy) 

Distance from 
cO (mas) 

Pos. Angle 
(°) 

FWHM 

(mas) 

2012.07 

55953.50 

cO 

708 ± 43 



0.08 ± 0.03 



d7 

34 ± 13 

0.13 ± 0.01 

-124.8 ± 11.2 

0.13 ± 0.04 



dlO 

346 ± 26 

0.13 ± 0.01 

-123.7 ± 8.6 

0.13 ± 0.03 



18 


2012.18 

55991.50 

cO 

756 

± 

45 







0.08 

it 

0.03 



d7 

66 

± 

16 

0.28 

it 

0.05 

-123.7 

it 

20.4 

0.16 

it 

0.03 



dlO 

264 

± 

22 

0.13 

it 

0.02 

-126.2 

it 

11.5 

0.14 

it 

0.03 

2012.25 

56019.50 

cO 

372 

± 

26 







0.06 

it 

0.03 



d7 

42 

± 

13 

0.41 

it 

0.08 

-122.8 

it 

14.2 

0.18 

it 

0.03 



dlO 

186 

± 

18 

0.13 

it 

0.01 

-113.1 

it 

6.1 

0.08 

it 

0.03 

2012.40 

56073.50 

cO 

680 

± 

41 







0.05 

it 

0.03 



d7 

33 

± 

11 

0.46 

it 

0.03 

-119.6 

it 

4.4 

0.07 

it 

0.03 



dlO 

195 

± 

17 

0.18 

it 

0.01 
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dlO 
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-112.7 
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dlO 
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0.03 



dll 
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0.46 

; it 

0.3 

-111.7 
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dlO 
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dll 
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d9 

dlO 

dll 

180 ± 19 
611 ± 32 

58 ± 14 

0.33 ± 0.03 
0.12 ± 0.01 

1.20 ± 0.07 

-107.9 ± 5.2 
-102.4 ± 3.2 
-114.6 ± 3.1 

0.18 ± 0.03 
0.10 ± 0.03 
0.20 ± 0.03 

2013.88 

56614.50 

cO 

788 ± 46 



0.05 ± 0.03 



d9 

367 ± 26 

0.26 ± 0.01 

-104.9 ± 1.3 

0.07 ± 0.03 



dlO 

677 ± 41 

0.12 ± 0.01 

-105.3 ± 1.6 

0.06 ± 0.03 



dll 

51 ± 15 

1.69 ± 0.10 

-116.8 ± 3.3 

0.26 ± 0.04 



dlla 

12 ± 1 

2.07 ± 0.14 

-114.9 ± 3.1 

0.17 ± 0.03 



dl2 

346 ± 26 

0.45 ± 0.01 

-109.5 ± 1.6 

0.13 ± 0.03 

2013.96 

56642.50 

cO 

623 ± 39 



0.09 ± 0.03 



dlO 

409 ± 29 

0.16 ± 0.01 

-113.9 ± 5.7 

0.14 ± 0.03 



dll 

58 ± 15 

1.72 ± 0.13 

-116.6 ± 4.0 

0.32 ± 0.03 



dlla 

36 ± 14 

2.22 ± 0.12 

-115.6 ± 3.7 

0.25 ± 0.04 



dl2 

223 ± 21 

0.59 ± 0.03 

-111.1 ± 2.6 

0.19 ± 0.03 



dl2a 

202 ± 19 

0.40 ± 0.02 

-109.6 ± 2.5 

0.13 ± 0.03 

2014.05 

56677.50 

cO 

639 ± 39 



0.07 ± 0.03 



d9 

98 ± 14 

0.27 ± 0.03 

-106.2 ± 5.2 

0.11 ± 0.03 



dlO 

322 ± 24 

0.13 ± 0.01 

-103.3 ± 4.1 

0.09 ± 0.03 



dll 

54 ± 15 

1.86 ± 0.15 

-115.0 ± 3.7 

0.36 ± 0.03 



dlla 

32 ± 14 

2.37 ± 0.15 

-116.0 ± 4.4 

0.29 ± 0.05 



dl2 

119 ± 17 

0.51 ± 0.06 

-108.5 ± 2.8 

0.22 ± 0.03 



dl2a 

79 ± 15 

0.51 ± 0.06 

-108.5 ± 4.2 

0.22 ± 0.03 

2014.15 

56713.50 

cO 

468 ± 30 



0.06 ± 0.03 



d8 

95 ± 14 

0.39 ± 0.03 

-108.6 ± 4.1 

0.12 ± 0.03 



d9 

109 ± 15 

0.23 ± 0.03 

-105.2 ± 5.7 

0.12 ± 0.03 



dlO 

338 ± 25 

0.08 ± 0.01 

-107.9 ± 4.9 

0.07 ± 0.03 



dll 

57 ± 16 

2.09 ± 0.19 

-114.5 ± 3.3 

0.46 ± 0.06 



dlla 

8 ± 2 

2 ± 0.12 

-115.1 ± 3.1 

0.12 ± 0.03 



dl2 

47 ± 14 

0.65 ± 0.13 

-111.8 ± 2.9 

0.20 ± 0.03 



dl2a 

20 ± 13 

0.65 ± 0.13 

-111.8 ± 4.2 

0.20 ± 0.03 

2014.33 

56780.50 

cO 

633 ± 38 



0.04 ± 0.03 



d8 

76 ± 14 

0.40 ± 0.03 

-110.5 ± 4.7 

0.12 ± 0.03 



d9 

146 ± 16 

0.24 ± 0.02 

-113.5 ± 4.3 

0.09 ± 0.03 



dlO 

463 ± 30 

0.10 ± 0.01 

-113.7 ± 2.3 

0.05 ± 0.03 



dll 

48 ± 15 

2.47 ± 0.12 

-114.2 ± 2.6 

0.29 ± 0.04 



dl2 

29 ± 13 

1.49 ± 0.09 

-111.4 ± 2.8 

0.18 ± 0.04 



dl2a 

54 ± 14 

0.87 ± 0.08 

-110.5 ± 4.8 

0.22 ± 0.04 



